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Abstract 

We present arguments for the existence of higher dimensional asymptotically AdS non Abelian black 
holes with a Ricci flat event horizon and analyze their basic properties. Unlike higher dimensional 
black holes with a curved horizon, of the usual Einstein- Yang-Mills system, these solutions have finite 
mass-energy. Below some non-zero critical temperature, they are thermodynamically preferred over the 
Abelian configurations. 

1 Introduction 

In recent years it became clear that various well-known, and rather intuitive, features of self-gravitating 
solutions with Maxwell fields in d = 3 -t- 1 spacetime dimensions are not shared by their counterparts with 
non Abehan gauge fields. For example, in contrast to the Abelian situation, self-gravitating Yang-Mills (YM) 
fields can form particle-like configurations [1]. The Einstein- Yang-Mills (EYM) equations also admit black 
hole solutions that are not uniquely characterised by their mass, angular momentum and YM charges [2]. 
Therefore the uniqueness theorem for electrovacuum black hole spacetimes ceases to exist for EYM systems. 
As a result, the literature on gravitating solutions with non Abelian fields has steadily grown up in the last 
two decades, including solutions with a cosmological constant A (see e.g. [3J|1] and references therein). The 
asymptotically anti-de Sitter (AAdS) solutions are of particular interest, since gauged supergravity theories 
playing an important role in AdS/CFT, generically contain non Abelian matter fields in the bulk, although 
to date mainly Abelian truncations are considered in the literature. Notably, non Abelian AAdS solutions 
exhibit new features which are absent for A = 0. For example, stable solutions with a global magnetic charge 
are known to exist even in the absence of a Higgs field [5] , [5] . 

In the context of the AdS/CFT correspondence [7], Klebanov and Witten have proposed a mechanism 
of spontaneously breaking gauge symmetry [8]. This mechanism has recently been exploited by Gubser et. 
al. to explain important phenomena in condensed matter physics, in particular superconductivity and 

critical phenomena. This mechanism results, with no recourse to supersymmetry, in a symmetry breaking 
boundary theory of a bulk gravitational theory with negative cosmological constant, the temperature of the 
black hole being nonzero. 

From our point of view, the most interesting development in this domain is the recent discovery in [9] that 
some AAdS non Abelian black hole solutions with a Ricci flat event horizon may posses superconducting 
horizons which are thermodynamically preferred below some non-zero critical temperature. Such solutions 
exhibit hair of the 'electric' part of the gauge field on the AdS boundary, manifesting the gauge symmetry 
breaking mechanism; at the same time the condensate of the 'magnetic' part floats above the horizon of the 
black hole. This mechanism was further exploited in subsequent works [10[ 111] . 

The only case discussed so far pertain to four dimensional AAdS spacetimes with Ricci flat horizon, and 
relatively little is known about such higher dimensional solutions with non Abelian matter fields. Naturally, 
it is always of interest to see how the dimensionality of the spacetime affects the physical consequences of 
a given theory. In particular, it would be interesting to see how general is the mechanism discovered in [9j. 
Besides, it is known from the work of JL3j, [14], [15], [16], that static spherically symmetric solutions of the 
usual gravitating YM system in spacetime dimensions d > 4 do not have finite energy as a result of their 
scaling propertes. Finite energy solutions exist only when the usual YM system is augmented with higher 
derivative corrections in the non Abelian action [16j , |17j . Therefore the examination of higher dimensional 



gravitating non Abelian solutions with a different topology (in this case with Ricci flat) of the event horizon 
is a pertinent task. 

Our objective in the present work is to extend this type of symmetry breaking to non Abelian EYM 
solutions used by Gubser in [9j, to arbitray d = D + 1 dimensional AAdS spacetimes. For this we exploit, 
qualitatively, previous results on d dimensional finite energy AAdS solutions for EYM systems given in [TB] . 
In fact, the actual model considered here differs from those of in that the latter is described by the purely 
'magnetic' components of the YM fielcQ whereas here our model will include also the 'electric' components 
which play an essential role as is the case of [5] . The salient feature of the YM models in j 16 is the presence 
of higher order terms in the YM curvature, whose role is to supply the necessary (Derick) scaling properties 
of the Lagrangian to enable the existence of finite energy solutions. There is however one major departure 
between the models in [16. and those exploited here. While the metric Ansatz employed in [TBI describes 
a spacetime with an S"'"^ x R boundary at infinity, the one here describes in contrast a flat Minkowkian 
boundary. One consequence of this is that the appropriate gauge group here is SO{D — 1), i.e. SO{d — 2), 
differing from the choice of SO{D), i.e. SO{d — 1), in [TB]. This feature is a reminder of the fact that the 
electric component Aq of the YM connection takes the role of a Higgs field. The other consequence of a fiat 
Minkowkian boundary is that inclusion of higher order YM curvature terms are no longer necessary for the 
solutions to describe finite energy configurations, as was the case when the black hole horizon had a nonzero 
Ricci tensor. This is a result of the much wider range of scaling properties satisfied when the metric has a 
Ricci fiat event horizon instead of the more restrictive scaling properties of the system when the horizon is 
spherical. Inclusion of higher order YM curvature terms, while not necessary for achieving finite energy, is 
still possible here, resulting only in quantitative effects. We have eschewed this option here since it is not 
qualitatively important. In addition, although we have restricted our attention here to spacetime dimensions 
5 < c? < 8 for simplicity, it is obvious that this limitation is unimportant. 

The metric Ansatz we use is a direct extension of that in [5J [HI HO] j to dimensions with a larger number of 
spacelike coordinates. To implement our procedure it is necessary to devise an Ansatz for the YM connection, 
generalising that used in previous work on YM fields in AdS spacetime. Here, we have found two distinct 
Ansatze which we have verified to be consistent. These generalise the distinct YM connection Ansatze of [9] 
and of [To], respectively. 

2 General formalism 

2.1 The field equations and the abehan solution 

Instead of specializing to a particular supergravity model, we shall consider the pure EYM theory with 
negative cosmological constant in d > 4 spacetime dimensions 

S = l d'x^g {j^{R - 2A) - i , (1) 

where the cosmological constant is A = — (d — 2)(rf — l)/2£^. Although it seems that the model ([T]) is 
non-supersymmetric in itself (at leasid for d > 4) , it usually enters the gauged supergravities as the basic 
building block. Therefore one can expect the basic features of its solutions to be generic. 

Variation of the action ([T]) with respect to the metric g^" and the gauge field leads to the EYM 
equations 

- ^.g^.i? + Ag^, = SttGT^,, D^F^''' = Q, (2) 

where T^. is the YM stress-energy tensor T^, = F^^F^^g'^^ - {g^^F^^F^'^'^ , and D^ = d^+ ig [A^, ■] (with 
g the gauge coupling constant). 

Inclusion of the 'electric' components of the YM field can readily be made, as e.g. in the case of Euclidean signature in 

m- 

^The case d = A, with A/(167rG) = —3g-^ corresponds to a consistent truncation of A/" = 4 gauged supergravity and may be 
uplifted to d = 11 supergravity [21], [22). [3l. 
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We shall consider black hole solutions of the above equation with locally flat horizons, which asymptoti- 
cally approach a locally AdS spacetime with a boundary at infinity R'^^^ . The simplest such configuration 
with a nonzero gauge field is represented by the Reissner-Nordstrom-AdS (RNAdS) black hole with 

'^^d ~ ^2 , fjcid-z) + ?'^^7^ ~ (§r ^ + Q2(j_2-)"' ^2(d-3) )'^^^; ^ — {c- ~ ^:^:^)Tdt, (3) 

■p- 73=3"'" g-'(d-2) ,.2(ti-3) 

where Mq, c, q are constants, T is an element of the gauge group while is the line element of the d — 2 
euclidean space. The parameters Mq and q are proportional to the mass and electric charge of the solution. 
The black hole horison is located at r = r^, with rl/l"^ - 2Afo/r^"^ + (87rG'(d - 3)/g2(d _ 2))q^ / r^^'^^^^ = 0. 
When taking q = 0, the Schwarzschild-AdS (SAdS) black hole with a planar horizon is recovered. 



2.2 The Ansatz 

We are interested in non Abelian configurations whose magnetic gauge potential vanishes asymptotically, 
such that the abelian configuration ([3]) is approached in that limit. Moreover, we shall suppose that our con- 
figurations present a dependence only on a suitable radial coordinate r which is orthogonal to the boundary 
of the spacetime. 

The choice of the gauge group compatible with these assumptions (and the corresponding YM Ansatz) is 
quite flexible. For the d = A case, two differenlH non Abelian Ansatze have been proposed in the literature, 
both of them for a gauge group SU(2). The first Ansatz used in [9^ corresponds to a "circular polarisation" 
of the magnetic YM connection and leads to an isotropic energy momentum tensor for the components on a 
surface of constant (r, i), with t the time coordinate. This is not the case for the second YM Ansatz proposed 
in [lOj . where a particular direction in the subspace is choosen, leading to a more complicated metric 
Ansatz. 

A straightforward generalisation of the isotropic Ansatz in f9| is found for a gauge group SO{D) (with 
D < d), the Ansatz for the YM connection being stated by 

A{r) = u{r)nl^T,jmidt + w{r)dx'T,D, (4) 

in which Tab — (r^ . Toy) are the gamma matrices in D dimensions and the indices a, P run over the range 
a,(3 = 1, 2, ..{D — l)/2. In the above relation, are the coordinates parametrizing a surface of constant r, t. 
Also, the sets ("■Q,"^^) appearing in ([4]) form a complete and orthonormal basis of constant valued vectors 
of unit length, 

D-1 D-1 D-1 -^5^ -^5^ 

i—l i—1 i—1 a — 1 a — 1 

The last completeness condition is very important and means that D — 1 is even. Therefore, ioi d = D + 1, 
the YM Ansatz ([4]) is definecQ only for an even number of space dimensions d. Of course, it can also be 
used for an arbitrary spacetime dimension by adding n codimensions y'^, with the non Abelian potential 
identically zero in the subspace labeled by the extra-coordinates. Obviously, this Ansatz leads to an isotropic 
energy momentum tensor in both and y -directions, with T^i ^ T\ . 

X y 

However, it is possible to define a different non Abelian Ansatz (Ansatz II in what follows), valid for 
all even and odd dimensions d. In this case we employ only one (constant valued) set of orthonormal unit 
vectors n*, in terms of which 

A{r) = u{r)n'Ti^d-idt + w{r)n^rjidx'^ . (5) 



^This contrasts with the case of solutions with a R X S'^ ^ boundary at infinity, where the choice of Ansatz is unique, see 
e.g. [6]. 

*The construction of a consistent isotropic Ansatz for d = 2k + I necessitates the enlargement of the gauge group, resulting 
in a YM connection with more than one magnetic potential. For example for the most interesting d = 5 case, the minimal 
gauge group is SO (5), while the Ansatz would contain six non Abelian potentials. 
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While the components of the connection ([5]) take their values fully in the algebra of SO{d — 1), only d — 2 
components of the YM connection are effective since nM^ = 0. Thus we have one additional radius, but the 
fields depend only on the radial coordinate r, so one of the magnetic components is zero. This results also 
on an anisotropic energy mementum tensor of the YM field. 

For an arbitrary spacetime dimension, a metric form compatible with the above two YM ansatze is given 

by 

d—n—2 n 

dsl ^ A{r)dr^ + Fi{r) ^ dx'dx' + F2{r)J2dy'dy' - B{r)df , (6) 

1=1 k=l 

where we have found convenient to take 

= B{r)=N{ry{r), F,{r) = f{r), = r2(/(r))2("+2-'')/" (7) 

with 

the function m(r) being related to the local mass-energy density up to some d— dependent factor. For the 
case of the first Ansatz Q, the YM field is defined on a subspace labeled by the (r, t; a;*)-coordinates, where 
i — 1, . . . ,d — n — 2 (with n > 0). The second YM ansatz corresponds to an arbitrary d, with n — 1 (Here 
we assume without any loss of generality n* = <5^_2 ^^'^ write x'^~'^ — .) 



2.3 The equations of motion and asymptotic solutions 

Within these Ansatze, the EYM field equations reduce to a set of five ordinary differential equations which 
can be expressed in a unified form a^ (where one takes c = 1/2 for the Ansatz I and c = l/{d — i) for 
Ansatz II) 

2/' cr' N'\ , / u2 {d-n-i)w'^\ w 



f r a N J \Na^ r^p J N 

„ (d-2 a'\ , 1 uw^ 

_ d-n-2 fNuP {d - n - , ^ , .9. 

2(7 f d — n — 2 2 n/2 



^=~pr\ — ^ ^ +a (w + 



2n ' ' ^~ ■ N2^2. 

■' ~ (d-n-2yf [n^ ^^2j2]^ / " y^^^o^y 

where = ?>TiG/{g'^{d — 2)). For d = 4, n = Q, the above equations reduce to those derived in [S] 
albeit for a different metric Ansatfl The Abelian Reissner-Nordstrom solution ([3]) is found for m(r) = 
Mo - c{d - 3)q^a^/2r''-^, f{r) = a{r) = 1, w{r) = 0, u(r) = uq - q/r<^-^. 

Unfortunately, no exact non Abelian solutions of this system are yet known. However, one can analyse 
their properties by using a combination of analytical and numerical methods, which are sufficient for most 



^Here and in what follows, the relations for the case with no codimensions are found by formally setting / = 1, followed by 
the limit n — > 0. 

= 4 AAdS non Abelian black holes with a Ricci flat horizon were discussed previously in |20| . 
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purposes. The solutions have the following expansiorj^ at r — > near the event horizon, located at r/j: 

m{r) = + mi(r - r^) + 0{r - rhf , (j{r) = (Th + oiir - rt) + 0{r - Vhf , 

I{r) = h + fi{r - rh) + 0{r - Th)\ (10) 
w{r) = Wh + 'Wi{r - Th) + 0{r - th)^ , u{r) = ui{r - rh) + U2{r - Thf + 0{r - r^f . 

Here, following |9j we interpret Wh as a 'magnetic' condensate. All coefBcients in the above relation can be 
expressed in terms of the real constants au, Wh, ,fh, ui- One finds e.g. 

rf-'^a^^ .4 4., ON 2 4^ . 2a2(d-n-3)nr^-5^2^2 



. . 2a,wW , {d-n-2) 2a^rl^-'tulwl \ 

' " ^ ' + n ^1"''^'^ + iid-l)ri-2m,rlPra, ) ' ^^^^ 

{d-n-3)fwl ui ( d-2 



' fi 

Note also that the physical condition N'{rh) > implies the following condition on the boundary data 
2c/>?/a2 + id-n- 3)wt/i4fi) < {d - \)l{aH^). 

We are interested in solutions of the EYM equations approaching at infinity the Abelian RNAdS solution 
([5]) . This implies the following asymptotic expansion as r — > oo 

- - ^^^{J' + c,H^)^ + 0(1A-^), air) = 1 - ^-^^ + 

f{r) = 1 - + Oil/r'^"-'^), w{r) = ^ + Oil/r"-'), u{r) ^ u, - + Oil/r'^"-'^), (12) 

with MQ,J,q,f real constants. The holographic interpretations of uo,q, and J are as follows: uq is the 
chemical potential, q is the electric charge, and J is that component of the current Ji on the boundary 
connected with the spontaneously broken part of the bulk gauge symmetry [9]. In other words, we have a 
D-dimensional conformal field theory described on the boundary of AdSo+i space equipped with SO{D — 1) 
conserved currents, which satisfy their own current algebra^. The normalisable boundary value of the 
'magnetic' field w{r) = J/r"^'^ + 0{l/r'^~^) corresponds to the vacuum expectation value of the boundary 
currents proportional to J arising after symmetry breaking, and, the existence of the horizon 'magnetic' 
condensate Wh. 

The case d = 4, n = considered in [9J is special from the point of view of the asymptotic expansion, 
since the finite energy requitements are compatible with a nonvanishing value of the magnetic potential at 
infinity. The condition w{oo) = is imposed there by requiring w to make a finite contribution to the norm 
of the non Abelian potential However, one can easily see from the field equation Q that for c? > 4, 
w{oo) 7^ results in a divergent value of the mass function m{r), which gives further justification to the 
choice p^ . 

The constant uq in the asymptotic expansion (|12p corresponds to the electrostatic potential $ = uo/g, 
while q fixes the electric charge density Qe = (d — 3)q/g. Other quantities of interest are the mass-energy 
density M, Hawking temperature T and entropy density S, 

The constant J which enters the asymptotics of the magnetic non Abelian potential w{r) corresponds to an 
order parameter describing the deviation from the Abelian solution. 

^The case of Ansatz II for d = 4 is special, as the near horizon expansions of w{r) and /(r) are different in that case: 
w{r) = Wfi + W2{r — Tfi)"^ + 0{r — ri^)^ , /(r) = fh + f2ir — Tf^)'^ + 0{r — r;,)^. However, the horizon data is still determined 
by (Th,Wh,fh,ui. 

®The dual CFT global current Ji is defined through Tr AiJi, where Ai here is the asymptotic YM connection. 
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Figure 1: The profiles of typical solutions are shown as a function of radial coordinate for Ansatz I (a) and 
Ansatz II (b). 

3 Numerical solutions 

3.1 Scaling properties and general features 

We start by noticing that the equations Q are not affected by the transformation: 

r — > Ar, m ^ X^^^m, £ —> X£, u ^ Xu, a a/X (14) 

while w, a and / remain unchanged. Thus, in this way one can always take an arbitrary positive value for a. 
The usual choice is a = 1, which fixefl^ the EYM length scale L = ^^-kG j {g^{d - 2), while the mass scale 
is fixed by = {^ttG / {g"^ {d — 2)Y'^^^^^^ /G. All other quantities get multipfied with suitable factors of L. 
However, in what follows, to avoid cluttering our expressions with a complicated dependence of (G, d), we 
take a unit value for a and ignore the extra-factors of g and G in the expressions of various global quantities. 

For solutions with a spherical event horizon, the event horizon radius and the value of the magnetic 
potential on the horizon are independent parameters (see e.g. [5]). This is not the case for the solutions 
here, in which case one can always set — 1 without any loss of generality. This is a consequence of the 
following scaling symmetry of the system ([9]): 

r Xr, w Xw, u — > Am, m X'^^^m, (15) 

while (T, / and the cosmological constant remain unchangecf^. 

The system ^ presents in addition two more scaling symmetries associated with the functions a and 
/ {e.g. a Act, u Xu, t Xt etc.). In the numerical procedure these symmetries are used to set 
ct(oo) = f{oo) = 1 and thus to fix the horizon values of the functions a and /. Together with the other 
symmetries mentioned above, this leaves us with three numerically relevant parameters: Wh, ui and the AdS 
length scale £. Since equations ^ are invariant under the transformation w —w, only values of Wh > 
are considered. 

The equations ^ with boundary conditions implied in turns by pop , (|12p have been solved numerically, 
using a standard shooting method. As expected, the properties of the solutions obtained for the two distinct 
YM Ansatze (jlj, ([5|) are rather similar and thus we have preferred to present them together. For the first 
case, families of solutions have been constructed in a systematic way for d = 4, 6, 8 with n — 0, and d ^ 5,7 

^These are the units usually used in the literature on EYM solutions . Note also that II14I I is a generic property of the 
EYM system, shared by solutions with a different event horizon topologies. 

1° The global quantities scale as follows: M \''-^M, T ^ XT, S ^ X'^'^S, Qe ^ X'^-^Qe, $ ^ A*, J ^ X'^-^J. 
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with n — 1. Several configurations with = 6, n = 2 have been constructed as well. When choosing instead 
the YM Ansatz (P, we have constructed solutions in d = 4, 5 and 6 dimensions. For every considered value 
of £, we could find regular black hole solutions for only one interval < Wh < wf^. The value of wf^ increases 
as £ decreases, Wh — corresponding to the RNAdS solution ^ . 

In all these cases, we noticed a number of common features. The behaviour of solutions for generic 
initial data is such that w —^ wq ed large r (in which case the total mass-energy diverges), or else 
there is a singularity at finite r. Given {wh,i), solutions with the right asymptotic behaviour ()12|) exist 
only for a discrete set of values of ui. As in the well known case of the Bartnik-McKinnon solutions [Ij, 
the solutions here are also indexed by the node number of the magnetic potential z«(r). It turns out that 
the configurations with nodes represent excited states whose energy is always greater than the energy of the 
corresponding nodeless configurations, and are therefore ignored in what follows. 

For all solutions the functions m(r), cr(r) and u{r) always increase monotonically with growing r. How- 
ever, /(r) and u'(r) feature a more complicated behaviour. Tyical solutions are presented in Figure 1 for 
both Ansatze. For sufficiently small LUh, all field variables remain close to their values for the Abelian con- 
figuration with the same r^. Significant differences occur for large enough values of ujh and the effect of the 
non Abelian field on the geometry becomes more and more pronounced. 

3.2 Thermal properties and superconducting horizons 

For all cases we considered, finite energy solutions were found only for values of the Hawking temperature less 
than a critical temperature Tc. As in the d = A case in [9] this Tc is, within numerical error, the temperature 
at which the RNAdS solution admits a static linearised perturbation, with nonvanishing but infinitesimally 
smal0 w. Also, our numerical results indicate that goes to zero for some critical value of the AdS length 
scale t, but the corresponding solutions do not appear to have a singular behaviour there. 

These features are shown in Figure 2 for several values of d. For a given dimension the part of the 
parameter space above the curve corresponds to the unbroken phase, where only Abelian solutions exist. 

In Figure 3, we plot several quantities which are invariant under the transformation (fT5|) as a function 
of the ratio T /Tc- AF there is the difference in the free energy density, M — TS, between a non Abelian 
solution and the the RNAdS solution with the same T and Qe- As usual, AF < means that the non 
Abelian solution is thermodynamically favoured. 

In all cases there is a second order phase transition with simple critical exponents, from RNAdS solutions 
to solutions with normalisable non Abelian condensates. We have verified that for T^, < T < Tc, (with 

^'^Our numerical code usually provided good quality solutions for T > Tc/2. 
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Figure 3: J/Qe and /S.F/Q'f "^^ are plotted as a function of T/Tc for d = 5 (Ansatz I) and d — & 

(Ansatz II) non Abelian solutions. 



always around 0.7 Tc), the solutions satisfy the universal relation J/Qe = ji/2\/^ ^ T/Tc, where ji/2 depends 
on the model. 



4 Further remarks 



In this work we have presented arguments that the d — A picture discovered in J)j is generic for the higher 
dimensional case as well. Considering several values of c? > 5, we have found evidence for the existence 
of a second order phase transition with simple critical exponents, from the (electrically charged) RNAdS 
solutions with a flat event horizon to non Abelian configurations with a nontrivial magnetic field. 

One should also note that the existence of these finite energy non Abelian solutions with d > 4 represents 
a surprise in itself, since it contradicts the expectation based on the no-go theorems in [13] , [14] , [15] , [E] . The 
negative results in the latter were proven for configurations with a spherical topology of the event horizon, in 
which case the electric potential necesarily vanishes when d > 5. Purely magnetic EYM solutions with finite 
mass were found by considering corrections to the YM Lagrangian consisting in higher order terms of the 
Yang-Mills hierarchy of the form Lp ~ Tr F{2p)'^ where F{2p) is the 2p-form p-fold totally antisymmetrised 
product of the SO{d) YM curvature 2-form F{2) (see e.g. [TT] for asymptotically flat configurations and 
[TO] . [TS] for solutions with a cosmological constant). Such systems occur in the low energy effective action of 
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Figure 4: Left: The profiles of a d = 4 dyonic non Abelian EYM solution with spherical horizon (continuous 
line) is plotted together with a RNAdS configuration with the same Hawking temperature and electric 
charge. Right: A number of quantities are presented as function of the magnetic charge for non Abelian 
monopoles and RNAdS solutions with the same temperature for c? = 4 solutions with £ = 1. We have ploted 
the Hawking temperature (multiplied with a factor of ten, for better visualisation), the ratio between the 
entropies of the Abelian and the SU(2) solutions, the black hole masses and the difference between the free 
energies. 



string theory and are, in some sense, the non Abelian counterparts of the Lovelock gravitational hierarchvl^. 
By contrast, we have verified that to construct finite energy EYM solutions with Ricci flat horizons, inclusion 
of higher derivative terms in the YM curvature to the action ([T]) is not necessary. For this reason we have 
restricted our attention to the usual gravitating YM system, as inclusion of higher order YM terms results 
only in (unimportant) quantitative changes. 

An important question is, whether it is only for solutions planar event horizon that the non Abelian 
solution is themodynamically favoured over the Abelian one? To answer this question, consider simply the 
case of the d — A dyonic SU(2) black holes with spherical event horizon topology originally discussed in Ref. 
[6]. These solutions are found within the Ansatz 

ds^ ^^ + r\de' +sm^ed4>^)^a\r)N{r)dt\ with 7V(r) = 1 - + J, (16) 

N(r) r £^ 

A = - ^u{r)T^dt + w{r)Tidd + (cot Ot^ + w{r)T2) siwOd^p^, 

where Tq are the Pauli matrices. Without any loss of generality, by using the symmetry p4p . one can 
set 47rG' = g = 1. The problem reduces in this case to a system of four coupled ordinary differential 
equations. The properties of these solutions including the boundary conditions and the asymptotic expansion 
can be founcP^ e.g. in [6J. The generic behaviour of the solutions is such that they have a nonvanishing 
magnetic charge Qm = 1 — w^(oo). Non Abelian solutions with Qm — 1 are found for special values of 
u'[rh) and have grr9tt 7^ The dyonic Abelian RNAdS solution with unit magnetic charge corresponds 
to w(r) = 0, u{r) = uq + q/r, a{r) = 1, m(r) = Mq — (1 + q'^)/2r. Considering again the case of a 
canonical ensemble, we have found numerical evidence for the existence of non Abelian solutions which are 
thermodynamically favoured over the Abelian ones. An example of such a situation is presented in Figure 4 
for an AdS length scale £ — 3. There, the Hawking temperature and the electric charge are T ~ 0.0053 and 



^■^ However, inclusion of members of this gravitational hierarchy turns out to be of no practical utility because subject to the 
symmetries imposed such terms with the requisite scaling properties vanish. 



^^The function p{r) in [6] corresponds to l/<T{r) in the Ansatz JTBJ. 



Qe ~ 0.751 for both solutions, while the masses are slightly different: M(RNAdS)~ 1.334 and M ~ 1.33 for 
the non Abelian counterpart. This implies AF < and thus the existence of a phase transitiorF^. 

One should also remark that, for c? = 4, the presence of a non Abelian electric field is not crucial for the 
existence of a phase transition between Abelian and non Abelian solutions. Setting u{r) = in the Ansatz 
(|16p . our numerical results indicate the existence of purely magnetic non Abelian configuration which are 
thcrmodynamically favoured over the abelian solutions with the same magnetic charge^ ^^nd temperature. 
In Figure 4b we plot a number of relevant qualities for a family of d = 4 EYM monopole black holes with 
rii = 1 and £ ^ 1 and the corresponding RNAdS solutions. The gauge potential w(r) is nodeless for all 
solutions there. The solution with Qm — corresponds to the Schwarzschild-AdS (SAdS) black hole with 
a spherical horizon. One can see that for A = —3 all non Abelian solutions with rh — 1 have AF < 0. 
However, the generic picture is more complicated, with a nontrivial dependence on £, r^- 

A study of these aspects is beyond the purposes of this work and will be presented elsewhere. 

We close by remarking that the asymptotic AdS structure of the spacetime is crucial for the existence of 
such solutions. As proven in [5S] for c? = 4, the asymptotically flat EYM solutions have no magnetic charge 
while their electric part vanishes identically. Moreover, by using the data in 2 , one can easily verify the 
difference between the free energy of a SU(2) hairy black hole and the Schwarzschild solution with the same 
temperature is always positive. 
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